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IDENTIFYING PRIMORDIAL SUBSTRUCTURE IN NGC 2264 

Paula S. Teixeira 1,2 , Charles J. Lada 1 , Erick T. Young 3 , Massimo Marengo 1 , August Muench 1 , James 
Muzerolle 3 , Nick Siegler 3 , George Rieke 3 , Lee Hartmann 4 , S. Thomas Megeath 1 , and Giovanni Fazio 1 

ABSTRACT 

We present new Spitzer Space Telescope observations of the young cluster NGC 2264. Observations 
at 24 /im with the Multiband Imaging Photometer has enabled us to identify the most highly cm- 
bedded and youngest objects in NGC 2264. This letter reports on one particular region of NGC 2264 
where bright 24 fim sources are spatially configured in curious linear structures with quasi-uniform 
separations. The majority of these sources (~60%) are found to be protostellar in nature with Class 
I sp ectral energy distr i butio ns. Comparison of their spatial distribution with sub-millimeter data 
fromlWolf-Cha se et al.l l|2003j) and millimeter data from lPeretto et all (J2005) shows a close correlation 
between the dust filaments and the linear spatial configurations of the protostars, indicating that star 
formation is occurring primarily within dense dusty filaments. Finally, the quasi-uniform separations 
of the protostars are found to be comparable in magnitude to the expected Jeans length suggesting 
thermal fragmentation of the dense filamentary material. 

Subject headings: infrared: stars — open clusters and associations: individual (NGC 2264) — stars: 
formation — stars: pre-main-sequence 



1. INTRODUCTION 

NGC 2264 i s an extended hie rarchically structured 
young cluster (jLada fc Lada 2003J associated with a gi- 
ant molecular cloud in the Monocero s OBI complex 
and l ocated at a distance of 800 pc (e.g. lDahm fe S imon 
2005). The cluster i s a v ery well studied region (e.g. 
iHerbid Il 954 lWalked ll954). displaying evidence of on- 
going sta r formation such as a plethora of molecular 
outflows ijMargulis et all 119881 IWolf-Chase et all 12003(1 
and Herbig-Haro objects l|Reipurth et al.l 12004ft . Lu- 
minous far infrared sources were detected by the In- 
fraRed Astronomical Satellite (IR AS), many of them be - 
ing Class and Class I sources ijMargulis et al.l ll989). 
However detailed study of the protostars in NGC 2264 
was hampered by the limiting sensitivity and resolution 
of IRAS. The Spitzer Space Telescope provides signifi- 
cantly higher resolution and sensitivity. We have thus 
re-visited NGC 2264 by surveying the cluster with both 
the Infrared Array Camera (IRAC) and the Multiband 
Imaging Photometer for Spitzer (MIPS). 

Here we present initial results of a MIPS and IRAC sur- 
vey of a very young star forming region wit hin NGC 2264, 
origin ally o bserved in the i nfrare d by iSargent et alJ 
l|1984ft and iSchwartz et"aTl l|1985fl . This particular 
region has been also obser ved at other wavele ngths, 
namely in the sub - millim eter (Willia ms fc Garl andl2002T: 
[W ^lf- Chase et aLl |2£)Q3) and millimeter l|Peretto et all 
2005). These observations reveal several filamen- 
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tary structures extending from a central IRAS source, 
IRAS 06382+0939 (a lso designated as IRAS 12 by 
iMargulis et alJ l|1989f) who classified it as a Class I 
source). IRAS 12's position appears to be a few arc- 
seconds away from a near-in frared binary discovered by 
iCastelaz fe Grasdalenl l|1988|) . R NO-West (B-typ e star) 
and RNO-East (low-mass star). S imon & Dahm (2005) 
have photometric and spectroscopic observations of this 
binary in X-ray and near-infrared wavelengths and pro- 
pose that RNO-E could be a proto-Herbig Ae star. 
Our Spitzer data indicate that IRAS 12 is in fact co- 
incident with this binary. Several sub-millimeter cores 
ijWilliams fc Garlandll200l IWolf-Chase i et al.ll2003^ and 
millimeter cores |peret^oet*^Ll ~ l2005fi are identified 
within these filaments. We compare results obtained 
with our Spitzer data with these sub- millimeter and mil- 
limeter observations in Sj3] Finally, we discuss primordial 
cluster sub-structuring in this region of NGC 2264 and 
comp are these results with those found for the Trapez- 
ium l|Lada et alJl2004|) in g| 

2. OBSERVATIONS AND DATA REDUCTION 

NGC 2264 was observed with IRAC as part of the 
Spitzer Guarant eed Time Observation program 37 
ijFazio et al.l2004|) . The data were acquired in two epochs 
seven months apart (March 6, 2004 and October 8, 2004), 
with two dithers at each epoch to allow easy removal of 
asteroids and other transients. The total mosaicked area 
corresponds to ~ 0.7° x 1.2°. The observations were 
performed in all IRAC bands (centered at 3.6, 4.5, 5.8 
and 8.0 /im) using the 12 s IRAC High Dynamic Range 
mode, consisting of two consecutive exposures with 0.4 
and 10.4 s integration time at each dither position. Ba- 
sic data reduction and calibration were done with the 
Spitzer Science Center pipeline, version S10.5. A final 
mosaic was created for each of the two HDR exposures, 
using the SSC mosaic software the MOPEX (version 
10/15/04), resampling the individual images on a final 
pixel scale of 0.86267"/pix (l/\/2 of the original IRAC 
pixel scale) to have optimal point source registering at 
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all IRAC bands. Cosmic rays and other outliers were 
removed using MOPEX temporal outlier module, and 
the diffuse background emission was matched between 
individual frames using the MOPEX overlap correction 
module. 

We generated a list of point sources in each mo- 
saic, at each wavelength , using the SExtractor package 
ijBertin fc Arnoutslll996|) . and then performed aperture 
photometry using the IRAF routine APPHOT with an 
aperture radius of 2 pixels (~ 1.73"). IRAC band mer- 
ging was performed using the IDL function srcor 5 , where 
we used a matching radius of 2 pixels. The saturated 
sources had their magnitudes replaced with those ob- 
tained by the short exposure ob servations (0.4s). 

The MIPS (jRieke et al.ll2004|) observations were con- 
ducted on 2004 March 16 using the scan map mode. 
Fourteen scan legs of 0.75° length and 160" offsets were 
taken at medium speed. Total integration times of 80 
s per point and 40 s per point were obtained in the 24 
/im band and 70 /im band, respectively. We also ob- 
tained sparse coverage in the 160 /im band, but most 
of those data were saturated due to the extremely high 
backgrounds from the molecular cloud. The full map was 
centered at 6 h 40 m 55 s +9°37'08" with a position angle of 
179°. These observations we re processed with th e MIPS 
Data Analysis Tool (DAT) {Gordon et all 12003) which 
produces calibrated mosaics of the mapped regions. Pro- 
cessing of the resultant image products to obtain pho- 
tometry was done using DAOPHOT and IDL routines. 
We used a zero point of 7.3 Jy for the 24 /im magnitude 
scale. 



IRAS 12, resembling "spokes" of a wheel. We therefore 
refer to this region as the "Spokes" cluster. We per- 
formed two simple statistical tests to ascertain whether 
the linear spatial alignments of sources could occur by 
chance. The first test consisted of examining 10,000 syn- 
thetic fields, with randomly generated distributions of 
objects. Each field has the same number of objects as 
the bright 24 /Km sources in the Spokes cluster. These 
fields were searched for linear alignments that consist of 
5 or more stars within a cone of an opening angle < 10° 
and a length < 3'. The second test consisted of radi- 
ally binning the bright 24 ^m sources, using IRAS 12 as 
the center, and comparing the average bin density with 
a Poissonian distribution. Both tests yield a probability 
of 0.01-0.02% of finding an alignment of stars like that 
observed in the Spokes cluster. 

To analyze the nature of the sources in the Spokes clus- 
ter we constructed Spectral Energy Distributions (SEDs) 
for all the sources in the selected field using Spitzer 
data. We determined the spectral indices, aiRAC, of 
those sources detected in all four IRAC bands (3.6 /im, 
4.5 /im, 5.8 /im, and 8/tm) by performing a linear fit to 
the IRAC flux points in a log(A) vs. log(AFA) diagram, 
for each source. Identifying Class I objects as sources 
with Chirac > (Lada 1987), we estimate the total frac- 
tion of Class I 24 /im sources in the entire field to be 33%. 
To better comprehend the spatial alignment of these 



3. RESULTS AND ANALYSIS 

Figure ^ shows a color composite image of the region 
in NGC2264 surrounding IRAS 12, which is the bright 
saturated source. The more extended nebulosity (green) 
in the image corresponds to polycyclic aromatic hydro- 
carbon (PAH) emission from the molecular cloud at 7.6- 
8 /tm and detected in the IRAC 8 /im band. The brightest 
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Fig. 1. — False color image of the Spokes cluster constructed 
from MIPS 24 ^m (red), IRAC 8.0 /im (green), and IRAC 3.6 /tm 
(blue) data. The image shows unusual linear spatial alignments 
of the brightest 24 /im sources. The central saturated source is 
IRAS 06382+0930 (IRAS 12). 

24 /im sources (red) form several linear structures that 
are not clearly evident in the shorter wavelength data. 
The linear arrangements of the brightest 24 /tm sources 
appear as arms extending from the brightest source, 

5 from IDL Astronomy Use r's Library: 

http://idlastro.gsfc.nasa.gov/homepage.html adapted from 
software from the Ultraviolet Imaging Telescope. 



Fig. 2. — Comparison of the spatial locations of 24 /im MIPS 
sources with dust emissi on at 850 fim (SCUBA data courtesy from 
IWolf-C hase et al. (200ij)). The greyscale and contours represent 
the sub-millimeter dust emission (contours range from to 2 in 
steps of 0.1 Jy/beam), while the star symbols mark the position of 
the sources detected at 24 /im with MIPS. The two squares mark 
the positions of the two brightest 24 fim sources (< 2 magnitudes). 
The sizes of the star and box symbols are proportional to the mag- 
nitude of the sources (ranging from 2.0 to 6.0 magnitudes for the 
star symbols). The beam size for the 850 fim data is indicated in 
the lower right corner. 



protostars we examine sub-millimeter observations by 
IWolf-Chase et all (|2003J) which trace dust emission from 
dense molecular gas. Figure 121 sh ows 850 /im e mission de- 
tected in the region IWolf-Chas e et al.l l2~003). The star 
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are spatially coincident with 9-10 of the sources listed in 
Table l 6 (other millimeter sources detected outside the 
filaments also have 1RAC and/or MIPS detections). 



Spectral index a IFiAC 

Fig. 3. — Diagram of the distribution of spectral indices 
for sources detected at 24 fim within (solid line) the 850 fim 0.4 
Jy/beam contour (see Figure and for the remaining sources in 
the field (dotted line). The vertical solid line corresponds to the 
median of the distribution of spectral indices of sources inside the 
contour, c«rRAC = 0-l, whereas the vertical dotted line corresponds 
to the median of the distribution of spectral indices of sources out- 
side the contour, ojirac = -1-1. 



symbols mark the positions of the bright 24 /im sources 
and the symbol sizes are proportional to the magnitudes 
of the sources. It is remarkable how well the bright 24 /mi 
sources are aligned with the dusty filaments. The same 
result is obtained when using the 450 /im data, also from 
IWolf-Chase et~aTl lE)0l . The comparison with the sub- 
millimeter data shows that the linear alignments have 
drawn our eye to filaments of the interstellar medium 
that appear to be breaking up and forming stars in a 
very regular pattern. 

We compare the distribution of Chirac f° r sources 
within the dusty filaments and for sources in the surroun- 
ding region, using the 850 /zm 0.4 Jy/beam contour in 
Figure El as the boundary separating these two regions. 
Figure [21 shows a histogram of the spectral indices de- 
termined for both the 24 /im sources within (solid line) 
and outside (dotted line) this contour. The diagram in 
Figure 01 shows that there is a higher concentration of 
Class I sources within the dusty filaments, as seen by 
comparing the median values of Chirac for both distri- 
butions: 0.1 for sources inside the filaments and -1.1 
for the remaining surrounding sources. The fraction of 
Class I 24 /im sources within the dusty filaments is 59% 
whereas for outside the filaments it is 23%. Not all the 
24 /im sources in the Spokes cluster have detections in 
all the IRAC bands. Of these remaining sources that 
have at least 3 band (IRAC and MIPS) detections we 
find 7 with rising slopes. We list all the Class I sources 
that we identified within the dusty filaments in Table 
1, sorted by their 24 /im magnitude. The spectral index 
Chirac is tabulated in column 4. We note that there are 
two 24 /im sources within the dusty filaments that have 
no IRAC counterpart, and as such we were not able to 
classify them. These sources could also be very young 
protostars. 

Recent millimeter observations bv lPeretto et alJ l)2005f) 
identify 1.2 mm emission peaks in this region that we find 



Neorest neighbor distonce (orcseconds) 

Fig. 4. — Histogram of the nearest neighbour distances between 
the 24 (im sources in the Spokes cluster (solid line) . We have re- 
stricted this sample to sources brighter than the 6th magnitude 
at 24 fim. The dashed line histogram corresponds to the averaged 
nearest neighbour distance distributions of 10,000 simulated fields 
populated with randomly positioned stars (see text for details). 
The preferential separation between the observed 24 /an sources is 
20" +5". 



Finally, we analyse the spatial separations of the proto- 
stars in the Spokes cluster. Figure 0] shows a histogram 
(solid line) of the nearest (projected) neighbor separa- 
tions between the brightest 24 fim sources in the Spokes 
cluster. There is a clear peak in the distribution corre- 
sponding to a separation of 20". The mean separation of 
sources in this distribution is 26". By using Monte Carlo 
techniques, we generated 10,000 synthetic star fields of 
the same area and number of stars as in the observed 
field. The dashed-line histogram in Figure 01 corresponds 
to the averaged nearest neighbor distributions of these 
simulated fields. We calculated the probability of gener- 
ating a field with the same fraction of sources, i.e. 44.7%, 
with nearest neighbor separations between 15" and 25" 
as in the Spokes cluster field. The probability of find- 
ing a field with a 40% or greater fraction of sources with 
those separations is 0.08%. The 20" peak is therefore 
v ery u nlikely a random featu re. 

IWilliams fc Garlandl l|2002|) determined a mean density 
from their sub-millimeter observations of 3xl0 4 cm~ 3 
for the dense gas in the Spokes cluster. Assum ing a 
temperature of 17 K ( Ward -Thompson et alJ I2000J) the 
corresponding Jeans length for this density is 0.105pc 
or 27" at a distance of 800 pc. For a mean density 
of 4.2xl0 4 cm -3 the corresponding Jeans length would 
be 20". These Jeans lengths agree very well with the 
separations we measure, possibly indicating that the 
dense filaments have thermally fragmentated into quasi- 
equidistant cores where the stars are now forming. This 

6 specifically D-MM1, D-MM2, D-MM3, D-M M6, D-MM8, 
D-MM 9, D-MM10, D-MM13, and D-MM15 from IPeretto et ail 
12005), corresponding to sources with # 4, 9, 7, 16, 5, 10, 15, 
11, and 2, respectively. 
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TABLE 1 

Class I sources in the Spokes cluster dusty 

FILAMENTS 



ID 


position (J2000) 
R.A. DEC. 


Q IRAC 


m [24 /im] 


1 


06 h 41 m 02.8 s 


+09° 36' 16" 




0.44 


2 


06 h 41 m 04.5 s 


+09°36'21" 


2.33 a 


1.00 


3 


06 h 41 m 05.8 s 


+09°35'29" 


0.52 


2.20 


4 


06 h 41 m 05.6 s 


+09°34'08" 


0.03 b 


2.21 


5 


06 h 41 m 01.8 s 


+09°34'34" 


0.13 


2.93 


6 


06 h 41 m 09.9 s 


+09°35'41" 


2.36 


2.94 


7 


06 h 41 m 04.3 s 


+09°34'59" 


0.09 


3.07 


8 


06 h 41 m 06.3 s 


+09° 33' 50" 


1.36 


3.33 


9 


06 h 41 m 06.5 s 


+09°35'54" 


3.26 b 


3.42 


10 


06 h 41 m 06.8 s 


+09°33'35" 


1.42 c 


3.54 


11 


06 h 41 m 08.6 s 


+09°35'42" 


1.99 b 


3.87 


12 


06 h 40 m 58.0 s 


+09°36'40" 


1.15 


4.11 


13 


06 h 40 m 59.1 s 


+09° 36' 15" 


0.65 b 


4.11 


14 


06 h 40 m 58.6 s 


+09°36'39" 


1.23 


4.19 


15 


06 h 41 m 07.7 s 


+09° 34' 19" 


2.22 


4.29 


16 


06 h 40 m 58.0 s 


+09°36'15" 


0.12 


6.77 


17 


06 h 41 m 04.5 s 


+09°33'44" 


0.02 


7.22 



a aiRAC determined by fitting 3.6 (im, 4.5 (im, and 8.0 fim 
b aiRAC determined by fitting 4.5 (im, 5.8 (im, and 8.0 /tm 
c c*lRAC determined by fitting 5.8 /im, 8.0 (im, and 24 fim 



res ult is reminiscent of that found in the Taurus clouds 
bv lHartmaiml <|2002j) . 

4. DISCUSSION AND CONCLUSIONS 

Our analysis of the Spitzer observations of NGC 2264 
leads us to conclude that we are identifying the primor- 
dial substructure of this cluster. Bright 24 /im sources, 
spatially arranged in unusual linear patterns, were found 
to be lying along and within dense fingers of molecular 
material with a spacing similar to that expected for sim- 
ple Jeans fragmentation. These sources are likely very 
young if we consider the short crossing times of the asso- 
ciated sub- millimeter and m illimeter cores suggested by 
Williams fc Garland! <|2002D of 0.5 Myr and lPeretto et all 
(2005) of 0.63 Myr, respectively. In fact, we find the ma- 
jority (~60%) of these sources to be Class I sources and 
protostellar candidates, suggesting that star formation 
in the Spokes cluster is occurring primarily within dense 
filamentary molecular structures. 

Primordial filamentary substructuring of this kind has 
also been fo und in the Trapezium cluster, in Orion. 
lLada et all (J2004), by analyzing deep ground-based 
3.4 /im observations of the Trapezium, find a deeply em- 
bedded population of young objects that traces a filamen- 
tary molecular ridge lying behind the main cluster. The 
older, more evolved foreground population in the Trapez- 
ium is distributed in a more dispersed manner with an 
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isothermal-like distribution. Similar to the Trapezium, 
NGC 2264 also exhibits a somewhat older and more dis- 
perse d population near the Sp okes cluster l)Lada fe; Ladal 
12003]) . iWilliams fc Garland! (12001 ob served H 13 CO+ 
emission lines and iPeretto et alJ l|2005fl observed N 2 H+ 
lines toward the sub-millimeter and millimeter peaks in 
the Spokes cluster, which they used to measure the core 
velocity dispersion and infer the virial mass of the sys- 
tem. They find that the virial mass is less than the total 
gas and dust mass enveloping the cluster. This means 
that the cores within the dusty filaments are currently 
bound to the filaments. If the total stellar mass of the 
cluster is less than the total mass of the dense gas in 
the filaments then the stellar group is destined to be- 
come unbound if the surrounding gas and dust is rapidly 
dispersed. This situation suggests that after star for- 
mation occurs in filamentary molecular fragments of the 
cloud, the more recently formed stars eventually expand 
to merge with the older population as the cloud material 
is rapidly cleared through outflows and jets, leading to 
the less structured d istributions of the overall population 
(jBonnell et al-lEoO^) . 

Finally, we find that the bright 24 fim sources ap- 
pear to have a preferential separation between themselves 
of 20" t 5". This distance is in very good agreement 
with the Jeans length corresponding to the mean density 
of the IRS-2 region determined bv IWilliams & Garlandl 
( 2002). It appears we are observing the result of the 
thermal fragmentation of the dense filamentary material 
into quasi-equidistant star-forming cores. On the other 
hand , turbulent motions a ppear to be present in this re- 
gion: IPeretto et alJ l)2005fl find that the millimeter peaks 
have velocity dispersions greater than the thermal sound 
speed for molecular gas at 15 K. Perhaps the dusty fila- 
ments may have formed thr ough turbulent motions of the 
cloud (Bonnell et al. 2003), then as the turbulence de- 
cayed the filaments became thermalized and fragmented 
into star-forming cores. The newly formed stars would 
generate outflows that re-energize the surrounding re- 
gion, feeding turbulence back into the cores and filaments 
and eventually disrupting the cloud. The positions of the 
young protostars may be the only remaining indication 
that thermal pressure played an important role in the 
formation of these objects. 
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